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ABSTRACT 

We discover that the mass of dark matter particles moM is imprinted in phase-correlations of the 
cosmic density field more significantly than in the 2-point correlation. In particular, phase-correlations 
trace moM out to scales about five times larger than the 2-point correlation. This result relies on a new 
estimator l(r) of pure phase-information in Fourier space, which can be interpreted as a parameter- 
free and scale-invariant tracer of filament-like structure. Based on simulated density fields we show 
how moM can, in principle, be measured using -^(r), given a suitably reconstructed density field. 



1. INTRODUCTION 

The model of a fiat and nearly scale-free universe dom- 
inated by dark energy and cold dark matter (ACDM), 
passed stringent empirical tests of the new millennium. 
The six free mo del parameters were found simultane- 
ously consistent (jKomatsu et al.ll2Qlll ) with the temper- 
ature fi uctuations in the c osmic microwave background 
(CMB, E arson et"aIll2Qlll ) measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP, 'Bennett et al.l 
[IbOSa.b), the baryon acoustic oscillations (BAOs) in the 
late-time large-scale structure (LSS) derived from galaxy 
redshift surveys (S DSS: Percival et al. 2010: 2dFGRS: 
iPercival et aTll2QQ7l : WigdeZ: lBlake et aLllMlt ). and dis- 
tance measurements based on type la supernovae (SNe, 
iHicken et aIll2QQ9l : iKessler et al.ll2QQ9[ ). 

This phenomenal success of the ACDM cosmology 
contrasts with our ignorance regarding the nature of 
its dark constituents. Crucial properties of these con- 
stituents, such as the particle mass of dark matter, are 
covertly imprinted in the sub-cluster structure o f the LSS 
(|Smith fc Markovid I2Q11I : ISchneider efall [2QT2i) . Yet, 
this information is not readily accessible to measure- 
ments. For one thing, the actual LSS is not directly 
observable due to the invisibility of dark matter, redshift- 
space d istortions (iKaiser 19871 ) ^ and general relativistic 
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effects (iBonvin fc Durrerl 

lYoo et anT2QQ9l ). For another, the information in the 
LSS is masked by a random component originating from 
quantum state reduction in the primordial universe. 

To filter out the random component, the observed 
LSS is usually subjected to statistical measures that 
are independent of cosmic randomness up to a volume- 
dependent shot noise, known as cosmic variance. In 
the case of a statistically homogeneous and isotropic 
universe, the infinite family of isotropic n-point cor- 
relation functions (n-PCFs) removes all randomness, 
but preserves all information (jFrvl H985). However, 
so far no finite set of statistical measures is known, 
which exclusively and exhaustively describes the infor- 
mation imprinted in the cosmic density field. Most cur- 



rent studies bypass this issue by considering only the 
isotropic 2-PCF ^2(^) or, equivalently, the power spec- 
trum p{k)^ where r and k denote the separation scale 
and wave- number. In doing so, important information 
is lost; e.g., subtle structural features, such as cosmic 
filaments, become indistinguishable from spherical fea- 
tures. Some studies improve on those drawbacks by 
invoking highe r-order co r relations (|Frv fc PeeblesI 119781 : 
ISuto fc Matsubaxal 119941 : iTakada fc Jain 1003) and al- 
ternative statistical measures, such as the fractal cor- 
relation dimension fScrimgeou r et al.l l2012f ). void dis- 
tribution functions (White 1979| ), and various shape- 
finders (Babul fc Starkman 1992; Luo fc Vishniac 1 9951: 
Sahn i et al. 1998; Aragon-Calvo et al. 2007; Bond e t al.l 
I2OIOI : ISousbiell201lf ). However, the benefit of these mea- 
sures in addition to ^2(^) is often limited, since they are 
heavily correlated to ^2(^) in terms of ensembles. To 
truly avoid this issue one must refer to statistical esti- 
mators that only measure information not yet contained 
in ^2(r) (e.g. Watt s et al.l 120031) . 

The aim of this work is to introduce a new statisti- 
cal estimator of the cosmic density field, which is based 
solely on the phases of the Fourier spectrum of the den- 
sity field, but not on its amplitudes, since the latter are 
already fully captured via ^2(^)- This requirement com- 
bined with the requirement of statistical homogeneity 
and isotropy naturally leads to a measure, which we will 
call the line-correlation function l{r). We show that, in a 
limited sense, this function can be interpreted as a proxy 
for cosmic 'filamentary' on length scales 2r. 

Unlike the 2-PCF, phase-correlations are independent 
of linear growth of LSS. Measures of phase-correlations 
such as l{r) are therefore particularly sensitive to the 
gravitational non-linear growth of the dark matter dom- 
inated cosmic web (, Watts et al.l 1200^ . In this work, we 
therefore chose to explore the dependence of i{r) on dif- 
ferent 'temperatures' of dark matter. Based on a series 
of large numerical A/'-body simulations, both with CDM 
and warm dark matter (WDM), we find that l{r) de- 
pends sensitively on the mass of dark matter particles 
moM- Our results suggest that ^{r) constrains moM an 
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order of magnitude better than ^2(^)- Moreover, £{r) 
depends on tubm on scales about five times larger than 
^2(^) - a pivotal result, since the complex baryon physics 
masking the footprint of dark matter properties becomes 
less important with increasing scales. 

The article proceeds as follows. Section [2] summarizes 
established concepts regarding cosmic structure and clar- 
ifies the meaning of n-PCFs and poly-spectra. Section [3] 
motivates and formally defines the line-correlation func- 
tion and presents geometrical interpretations. A range 
of cosmological applications, namely the measurement of 
mDM, is then considered in Section HI based on a series of 
A^-body dark matter simulations. Section [5] summarizes 
the key results and discusses their potential application 
to observed data. 

2. COSMIC STRUCTURE AND CORRELATION 
FUNCTIONS 

This section reviews the statistical nature of cos- 
mic LSS and summarizes the concepts of correlation 
functions and spe ctral analysis (details in section 3 of 
iBernardeau et al.l [20020 . 

2.1. Cosmic density field and its statistical symmetry 

We consider a flat Euchdean universe, consistent with 
BAO measurements assuming a cosmological constant 
(jPercival et al]l2010l : rBlake et aLll201l[ ). with a mass den- 
sity field 



(1) 



where m denotes the mass, V the comoving volume, and 
r G the position in D spatial dimensions; for illus- 
trative purposes we consider both D = 2 and D = 3. 
To simplify the notation, we omit the implicit time- 
dependance of p in eq. ([1]). 

According to the Big Bang theory, p{r) evolved from 
a dense, maximally symmetric state under the action 
of physical laws that are spatially homogeneous and 
isotropic. Complex substructure then grew from seeds 
of reduced symmetry, known as quantum- fluctuations, 
caused by quant um state reductio n in the inflating pri- 
mordial universe (jLeon et al.ll201lf ). In the current view, 
quantum state reduction decreases the spatial symme- 
try, but maintains homogeneity and isotropy in the 
sense that the outcome probabilities of the proces s con- 
serve the symmetry of the evo lution operator (Mc Weenyl 
120021 : lObreschkow et al.ll2007[ ). This weaker, probabilis- 
tic symmetry is referred to as statistical homogeneity and 
isotropy. The resulting conjecture that p(r) is statisti- 
cally homogeneous and isotropic is known as the 'cos- 
mological principle' and is supported b y modern redshift 
surveys (e.g. SDSS: Gong 2010; Svlos Labini & Barvshev 
120101 : WiggleZ: Scrimgeour et al. 2012). We express the 
statistical homogeneity and isotropy explicitly by writing 
p(r) as 
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where G are random translation vectors and 

Ri G 0{D) are rotation matrices of the orthogonal group 
{det{Ri) = ±1). The generating functions ^i(r) > are 
defined such that all cosmological information is encoded 
in ^i(r), while all quantum randomness is absorbed in 



the variables ti and Ri. By definition, eq. (|2]) thus sep- 
arates non-random variables {^^(r)} from random ones 
{ti, Ri}. This separation is useful when constructing sta- 
tistical measures that isolate the information. 

Given the compelling observational evidence for the 
large-scale homo geneity, thus non -fractal structure, of 
the universe (e.g. lScrimgeour et al.l l2012) we can define a 
universal average density p and the density perturbation 
field 



p{r) - p 



> -1. 
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This field then satisfies (5y ^ as V ^ oo, with dy = 
J dV5{r) being the average density perturbation. 
Let us comment on a few points. First, a corol- 
lary of the spatial homogeneity is that the total mass 
fy dVp{r) = Vpy is proportional to F as V ^ oo. In 
other words, the fractal dimen sion converges to 3 in this 
limit (|Scrimgeour et al.ll2012[ ). Second, it is a common 
misconception that statistical homogeneity and isotropy 
only concern large-scale (> lOOMpc) averages. For ex- 
ample, a universe with all mass concentrated around the 
nodes of a Cartesian grid with 1 Mpc spacings would 
satisfy py ^ p as V ^ oo, but violate statistical ho- 
mogeneity as stated in eq. ([2j). Third, although our uni- 
verse seems statistically homogeneous and isotropic, ob- 
servational proxies of p(r), such as redshift-surveys, can 
violate this statistical symmetry. A famous example is 
the finger s-of- God effect ( Kaiserl 119871 ) originating from 
a Doppler-shift contamination in the observed redshifts 
- a potentially us eful feature for ob s ervati onal cosmology 
as emphasized by iRaccanelli et al.l ()2012f ) . 

2.2. n-point correlation functions 

Given a density field p(r) [eq. (|2])] that mixes informa- 
tion with random translations (homogeneity) and rota- 
tions (isotropy), how can we extract the information? 

In a first step, statistical homogeneity is exploited by 
averaging over all translations. This is the key idea be- 
hind the correlation functions, which are spatial averages 
of product functions. The n-point density correlation 
function (n-PCF) is defined as 



S,(ri,...,r,_i) = - y d^t []j(t + r,). 



(4) 



where = 0. In particular, the 2-PCF reads 

S2(r)^i Jd^tSimt + r). (5) 

In a second step, statistical isotropy is exploited by av- 
eraging over all rotations R G 0{D). This leads to the 
'isotropic' n-PCFs, 



^n(5{ri,...,rn-i}) = En{Rri,...,Rrn-i)ji, (6) 

where <S{ri, r^-i} denotes a unique representation of 
the shape defined by the n-point set {0,ri, ...,rn-i} re- 
gardless of its orientation. In the case of n = 2, this 
shape reduces to the distance r = |ri|. The resulting 
isotropic 2-PCF 



6(r)^S2(i?r), 



(7) 



is by far the most common statistical measure of LSS, as 
justified in Section [2^ 

We emphasize that and here refer to single 
realizations of the density field and not ensembles of 
fields, i.e. En ^ (S^) and 7^ (^n), where ( ) de- 
notes the ensemble average. Furthermore, these n- 
PCFs refer to the perturbation field 5{r) rather than 
p(r). Our 2-PCF and 3-PCF are therefore identical to 
those called the 'reduced' 2-PCF and 3-PCF by Peacoc2 
(1999) and t he 'connected parts' of 2-PCF and 3-PCF 
by Bernardea u et al.l (j2QQ2l ). The family of the isotropic 
n-PCFs is statistically complete ( Frvl Il985f ) in that 
it contains all the information contained in the density 
field p(r). This information is contaminated by the cos- 
mic variance ((£ n ~ (£ n))^), which can be calculated for 
any function (iSzap udi 2001) and vanishes as ^ ^ oo. 

2.3. Fourier space representations 

Since the correlation functions are convolution inte- 
grals over t G R^, they can be computed more efficiently 
in Fourier space. Using the standard Fourier transform 
FT : S{r) (5(k) in cosmology and its inverse (IFT), 

and expressing all (5(r) in eq. (j4]) as IFT[(5(k)], we find 
(details in Appendix R|) 



Sn(ri, ...,rn-i) 



n-1 



X J d^'kn-ie^'^-^-'-^ Pn(ki,...,kn_i), 



(8) 



where k G is the wavevector and the complex- valued 
functions 

Pn(ki, K-i) = (5(ki) • • • (5(kn_i)(5(-Sk,) (9) 

are called 'poly-spectra'. Thus, for any n > 2, the corre- 
lation function En is equal to the IFT (generalized to n—1 
variables) of the poly-spectrum P^. The most common 
poly-spectra are the (real) power spectrum P(k) = P2(k) 
and the bi-spectrum P(k, q) = P3(k, q), 

P(k) = ,5(k),5(-k) = |,5(k)|', (10) 

B(k,q) = ^(k)^(q),5(-k-q). (11) 

Like in the case of the n-PCFs, the spectra P(k) 7^ 
(P(k)) and i?(k, q) ^ (i?(k, q)) here refer to a single 
density field, not to ensembles thereof. These spectra 
are the FTs of the 2-PCF and 3-PCF, respectively. 



52(r): 



(2 



J 



k e'^ ' 



(12) 



^3(r, s) = 1 1 d^k d^q e'(k-r+q-s) B(k, q). (13) 

The relations between En and P^ imply similar rela- 
tions between the isotropic correlation functions and 
rotationally symmetrized poly-spectra, called 'isotropic' 
poly-spectra. 



Pn(^{ki,...,kn-l}) = Pn(i^ki,...,Pkn-l)^. (14) 

Fig. [T] depicts the hierarchy from the density field down 
to isotropic correlation functions and their equivalents 
in Fourier space. Because of their importance explicit 
expressions for ^2 and ^3 are given in Appendix [Bl 
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Fig. 1. — (Color online) Hierarchy of density field, correlations 
functions, and isotropic correlation functions together with their 
spectral equivalents. The irreversible mapping labeled 'homogene- 
ity' removes the random translations, which contain no informa- 
tion if 5{y) is statistically homogeneous. Similarly, the irreversible 
mapping labeled 'isotropy' removes the rotations, which contain no 
information if Sir) is statistically isotropic. 



2.4. Cosmological importance of ^2{r) 

In the c u rrent view (jBennett et al.l I2Q11I : 
iKomatsu et al.l I2Q11I ), the primordial density fluc- 
tuations imprinted in the CMB are consistent with a 
Gaussian random field (GRF). A GRF results from 
a random (Poissonian) superposition of infinitely 
many plane or spherical waves with vanishing phase- 
correlation. The evidence for the Gaussianity of the 
CMB thus supports the physical interpretation that the 
primordial density fluctuations derive from de-correlated 
quantum fluctuations producing a bath of spherical 
sound waves. A key property of a GRF is that its in- 
formation, i.e., its non-randomness, is entirely contained 
in the isotropic 2-PCF ^2(^)- Thus, as far as current 
measurements can tell, all the cosmological information 
of the CMB is contained in ^2(^), or, equivalently, in its 
isotropic power spectrum p{k). 

3. PHASE-INFORMATION AND 
LINE-CORRELATIONS 

This section introduces the 'line-correlation' function, 
a new estimator of phase-information of cosmic structure. 

3.1. What is phase- information? 

For a statistically isotropic density field S{r)^ ^2(^) in 
eq. ([7j) contains the same information as the full 2-PCF 
E2{r). The reversible mapping between E2{r) in eq. (fT2]) 

and the amplitudes |(5(k)| then implies that ^2(1*) mea- 
sures all the cosmological information contained in the 
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amplitude field |(5(k)|. All additional information, not 
captured by ^2(^)7 niust therefore reside in the phases- 
factors 

e(k) = ii^^=e^-s[^(k)l. (15) 
\5{k)\ 

The information contained in these phase-factors is 
called phase-information and it can take the form of 
phase-phase correlations (e(k)e(q)) and/or amplitude- 
phase correlations (|(5(k)|e(q)). Unlike the primor- 
dial universe, the local universe does indeed con- 
tain a significant amount of phase-information, as ev- 
idenced by its clearly non-vanishing isotropic 3-PCF 
(e.g. 2dFGRS: Croton et al. 2004; Gaztahaga et al. 2005; 
SDSS: Nichol et al. 2006; Marin 2011; numerical simu- 
lations: Barr iga fc Gaztaha ga 2002i) . Therefore, phase- 
information measurements of the la te-time LSS prom ise 
to be a pivotal cosmological probe (|Watts et al.ll2003l ). 



3.2. Concept of line- correlations 

A natural way to measure phase-information is to 
use higher-order correlations (n > 3). However, 
this choice may be problematic since the estimators 
(n > 3) and ^2 are strongly correlated in the sense that 
they correlate, i.e. cov(^2,^n) 7^ 0, even if there are 
no amplitude-phase correlations, i.e. (|(5(k)|e(q)) = 0, 
across the considered ensemble of density fields. This 
strong correlation between ^ 3) and ^2 is due to 

the fact that the poly-spectra Pn depend directly on the 
amplitudes |(5(k)| in addition to the phases e(k). Thus, 
constraints on cosmological parameters obtained from 
(n > 3) are generally strongly correlated to constraints 
obtained from ^2 , which can lead to serious statistical dif- 
ficulties. Alternative estimators, which are fully defined 
by phase-phase correlations (e(k)e(q)) without explicit 

dependence on the amplitudes |(5(k)|, are here called 
measures of pure phase-information. Such estimators, 
must be defined exclusively upon the phase-factors e(k), 
or, equivalently, on the field e(r) = IFT[e(k)]. 

Fig.[2]shows an example of a density field S{t) with the 
corresponding fields S{k) = FT[5(r)], e(k) = J(k)/|(5(k)|, 
and e(r) = IFT[e(k)]. In this example, ^(r) is a statis- 
tically homogeneous and isotropic 2D density field con- 
structed on the basis of eq. (|2]) . The generating functions 
g{T) are identical elongated 2D Gaussian distributions. 
The mapping 5(k) e(k) removes all information stored 

in the amplitudes |(5(k)| and therefore all 2-point corre- 
lations. Thus, the filamentary structure of e(r) shown in 
Fig. [2] exclusively represents phase-information of ^(r). 

By construction, e(r) exhibits vanishing 2-point cor- 
relations for all r 7^ an thus the simplest meaning- 
ful measure of pure phase-information of ^(r) must be 
based on 3-point correlations of e(r). As illustrated in 
Fig. [21 the removal of 2-point correlation tends to col- 
lapse elongated structures to fine segments. Therefore, 
the most natural 3-PCF to consider is that of three points 
on a straight line. For simplicity we chose these points 
to be equidistant. Using the explicit expressions for 
^3(r) = S3(r, — r)|^|^^ given in eq. (|B6[) and substitut- 



ing (5(k) for e(k), we then obtain the modified 3-PCF 

^3W=(^//d"^d"'^--(l'^-^l^)||^' 

with the kernel 

= isin(x)/x, ifD = 3. (1^) 
According to eq. ([8]), Q{r) is identical to 

e3*(r)=e(t)e(t + r)e(t-r),_|r|=,. (18) 

Unfortunately, ^|(r) is an ill-defined function. (Math- 
ematically, it is a distribution.) This can be seen when 
expressing e in eq. (fTSj) as a discrete IFT. To do so, we 
adopt the standard discretization scheme explained in 
Appendix [D] and illustrated in Fig. [TTJ In this scheme 
the density field is represented on a squared {D = 2) 
or cubic {D = 3) box with side- length L, A^^ grid cells 
spaced by Ar = L/N ^ and periodic boundary conditions. 
Its Fourier space becomes a regular box of side-length 
27r/Ar and cells spaced by AA; = 27r/L. The IFT 
then reads e(r) = ^(^)- ^ particular discretiza- 

tion is fully defined by L and Ar, and any physical func- 
tion /(r) should converge both as Ar and L ^ oo. 
Yet, both requirements are violated by Q{r). 

First, as Ar 0, the number of independent modes 
with wavelengths shorter than a fixed r diverges as 
Ar~^, since max{|k|} = ir/Ar. The ^ Ar~^ modes 
with wavelengths shorter than the characteristic scale 
of ^(r) exhibit vanishing amplitudes and thus random 
phase-factors e(k). Therefore, |e(r)| = | e*^"'^ e(k)| cx 
Ar"^/^, if Ar ^ . According to eq. ([18]), ^|(r) hence 
diverges as Ar~^^/^ and becomes infinitely dominated 
by random noise. To avoid this divergence, we must 
limit the number of high-frequency modes in eq. (fT6|) . By 
virtue of the Nyquist-Shannon theorem, it is natural to 
impose |k| < 7r/r and |q| < 7r/r. In signal processing ter- 
minology, this mode-truncation is called a low-pass filter 
with a spherical top-hat kernel. With this modification 
^l(r) becomes independent of Ar, if Ar < r/2. 

Second, as L oo, we face a similar challenge: the 
density of modes per volume of Fourier space increases 
as due to the mode spacing Ak = 27r/L. If L is 
longer than the longest physical correlation lengths, then 
increasing L corresponds to adding sub-modes with un- 
correlated phases. Hence, we are again in the random 
phase case, where |e(r)| = | e*^'"" e(k)| cx L^/^; ^|(r) 
then diverges as L^^/^. To avoid this, while keeping 
^3(r) dimensionless, we must divide ^Kr) by (L/r)^^/^. 

3.3. The isotropic line- correlation function i{r) 

Based on the conceptual discussion of Section 13. 2| we 
now define the 'isotropic line-correlation function' of the 
density perturbation field 5{r) as 

|k|,|q|<7r/r 

where wd{x) is given in eq. ([T7|) . Like the bispectrum 
5(k, q), given in eq. ([TT]) . the line-correlation function 
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Fig. 2. — (Color online) Key idea of Section [321 to find an estimator of the density field S(r) that is uncorrelated to ^2(^) up to residual 
correlations stemming from amplitude-phase correlations {|5(k)|e(q)}, we remove all 2-point correlation from S(r). In Fourier space this is 
equivalent to suppressing all amplitude information, i.e., to the mapping S(k.) e(k) = 5(k)/|5(k)|. Any statistical measure depending 
only on e(k) or e(r) = IFT(e(k)) is then de-correlated from ^2(t) for an ensemble of fields with vanishing amplitude-phase correlations. 
In this example, 6{r) is a superposition of 200 randomly shifted and rotated 2D Gaussian distributions eyip(—x'^/2a'^ — y^/2(7^) with 
(o'xTO-y) = (0.006L, 0.024L). The fields S(r), ^(k), e(k), e(r) are discretized using the scheme of Appendix [D] with N = 300 cells a side. 
The complex fields ^(k) and e(k) are represented with brightness for amplitudes and hue-colors for phases. 



£{r) ^ {£{r)) here refers to a particular realization of 
6(r) rather than an ensemble of fields. If some values 
of (5(k) vanish, then |5(k, q)| = and i{r) diverges. In 
practice, this only happens for k = 0, where <5(k) strictly 
vanishes by virtue of S{r) = 0; for this case we adopt 
e(0) = ^(0)/|^(0)| = 0. The discretized version of the 
function i{r) for the case of a periodic Cartesian grid 
with side- length L reads (see Appendix [D]) 

= ^ (20) 

|k|,|q|<7r/r ' ^ ' ' 

Given the definition of £{r) in eq. (fTOj) . a list of basic 
properties of can be derived: 

(i) £{r) measures statistically homogeneous and 
isotropic information in the field S{r). 

(ii) £{r) correlates to the 2-PCF and power spectrum 
only through amplitude-phase correlations, thus 
cov(^2,^) = if (|(5(k)|e(q)) = 0. In particular, 
for an ensemble of GRFs with random phases, ^2 
and i are statistically independent - a statement, 
which is not true for ^2 and > 3). 

(iii) £{r) is invariant with respect to addition S{t) 
6{t) + c and multiplication 6{t) cS{t) (c ^ 0), 
where c is a real constant. By virtue of eq. (j3j), £{r) 
of ^(r) is therefore identical to £{r) of /)(r). 



(iv) The mapping ^(r) £{r) is non-linear with respect 
to the superposition of two density fields. 

3.4. Physical interpretation of £{r) 

The correlation function £{r) is a measure of pure 
phase-information, i.e. it depends only on phase-phase 
correlations (e(k)e(q)). In turn, the 2-PCF ^2(^) depends 

only on amplitude- amplitude correlations (|(5(k)| |(5(q)|). 
The latter is often used as a measure of clustering. Sim- 
ilarly, £{r) can be interpreted as a measure of elongated 
structures, such as cosmic filaments, in a sense specified 
in the following. 

Fig. [3] shows four examples of a 2D random field ^(r), 
constructed by using eq. ([2j): (a) a GRF made of plane 
waves, (b) a random field with circular kernels, (c) a 
random field using filament-like kernels, and (d) a su- 
perposition of the GRF (a) and the filamentary field (c). 
The fields (a)-(c) are designed to exhibit the same 2-PCF 
^2(^)5 thus illustrating that ^2(^) cannot distinguish be- 
tween wave-like, sp herical, and filamentary substructure 
(see also Fig. 1 bv lColesll2QQ5f ). In contrast, ^3(r) van- 
ishes for the GRF, but still exhibits a similar shape for 
the spherical and filamentary fields. In principle, the full 
isotropic 3-PCF ^3(|r|, Zr, s) can distinguish spherical 
from filamentary structure (e.g. Nichol et al. 2006), but 
its geometrical interpretation is difficult. Only the line- 
correlation function £{r) clearly separates the filamentary 
density field (c) from the fields (a) and (b). Moreover, 
the shape of £{r) describes the straight filaments quanta- 
tively: £{r) exhibits a bump, £{r) > 0.2, roughly on the 
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Fig. 3. — (Color online) Example of three 2D density fields 6(r) with the corresponding correlation functions ^2(^)5 ^3(^)5 ci^d £(r). The 
shaded envelopes of the solid lines represent 67% confidence intervals. The density fields are chosen such that they all yield a similar 
^2(^) and such that two of them yield a similar ^3(r). (a) GRF constructed by superposing 350 randomly oriented plane waves; (b) a 
superposition of 200 symmetric 2D-Gaussian distributions exp(— r^/2cr^) with a = 0.014L; (c) a superposition of 30 randomly oriented 
asymmetric 2D Gaussian distributions exp(— x^/2crj — y^/2(jy) with ((Jx,cry) = (0.006L, 0.024L); (d) density field is equal to the sum of 
the density fields in examples (a) and (c). 
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interval r G [(Jxi(^y\ between the short and the long char- 
acteristic filament radius (see caption of Fig. [3]) . Hence, 
for a density field composed exclusively of straight fila- 
ments, l{r) measures the filamentarity on length scales 
2r. This feature of l{r) remains true, even if the filamen- 
tary field (c) is superposed with the GRF (a), such as 
illustrated in Fig. [3]i. 

A more systematic account of how aspherical substruc- 
ture is imprinted in ^{r) is provided in Fig. IH This anal- 
ysis relies on a 3D density field constructed via eq. (j2j), 
with the generating functions ^(r) being 50 identical 
spheroids. These spheroids are fully characterized by 
their prolateness defined as the ratio between the pole- 
radius Vp and the equator-radius rg , and by their average 
radius ro = (r^Tg)^/^. The left panel in Fig. |4] shows a 
projection of the 3D density field for the case of ^ = 8. 
The middle panel shows £(r) for three selected values of 
confirming that ^{r) is positive for r > min{re,rp} and 
peaks in the interval r G [re,rp]. For r < min{re,rp}, 
i{r) undergoes a series of oscillations. In the case of per- 
fectly spherical {q = 1) substructure i{r) does not vanish; 
however, the integral J i{r)dr nearly vanishes for g = 1 
and increases for both q < 1 and q > 1. We thus see 
that £{r) is particularly sensitive to aspherical substruc- 
ture on scales ~ 2r, be it oblate, such as cosmic sheets, 
or prolate, such as the more common cosmic filaments. 

The line-correlation function i{r) can also be used to 
characterize a mixture of aspherical structures. An il- 
lustration for the case of a 2D density field is shown in 
Fig.m This density field is again constructed via eq. (|2]), 
by randomly superposing two sizes of filaments. In this 
case i{r) is double-peaked with each peak approximately 
measuring the width and length of one type of filament. 

Next, we investigate the dependence of £{r) and 
J i{r)dr on the size and density of substructure, e.g. the 
length and density of filaments. Let us consider a density 
field constructed via eq. (|2]) by superposing m identical 
generating functions ^(r) of a characteristic linear scale 
ro and a characteristic volume Vb, e.g., spheroids with 
^0 = {rpviy^^ and Vb = (47r/3)rQ. We define the filling- 
factor of this density field as / = mVo/L^. If the filling 
factor is increased by increasing the number of objects 
m, then the increasing number of random translations 
and rotations in eq. (|2]) amplifies the phase-noise as ^/rn. 
Hence the amplitude of the line-correlation £{r) falls (on 
average) as 1/ ^/m. If, on the other hand, the character- 
istic scale ro is changed by a factor u while maintaining 
/ constant, i.e., by varying the number of objects per 
unit volume, the situation looks as follows. Let fo = u^o 
be the new length scale, and ^(r) the corresponding den- 
sity field. Since £{r) is independent of the box size L as 

L ^ oo, we may choose the new density field S{t) to be 
defined on a cubic volume with side- length L = uL^ such 
that the number of objects rh in the total volume remains 
the same, i.e., rh = m. In this way we find 5{ur) = S{r)^ 
given the same choice of random translations and rota- 
tions for 6{r) and S{r). The line-correlation associated 

with S then becomes £{r) = £{r/u). As a consequence, 
if in a density field S{r) with line-correlation £{r) the 
substructure-scale ro is stretched by a factor u and the 
filling factor / is varied by a factor the line-correlation 



scales as 

£{r) = v-^l^ £ (0 (21) 

It follows that the rescaled function \/J £{r/rQ) is inde- 
pendent of / and ro (up to shot noise) and thus character- 
istic of a particular type of substructure, e.g., spheroids 
of a certain prolateness q. This also implies that 

r-^ j £{r)dr = nf-^^^ (22) 

with n being a constant depending on the shape of the 
substructure, but not on its scale and filling factor. This 
relation is strongly supported by the numerical example 
in Fig. [6l 

In summary, £{r) is sensitive to aspherical substructure 
on scales 2r. In the special case of a filamentary field, 
£{r) measures the characteristic scales of the filaments. 
The dependance of £{r) on the prolateness of aspherical 
substructure is shown in Fig. HI The dependance of £{r) 
on the size and filling factor of substructure is given in 
eq. (|2T]) . and illustrated in Fig. [6l 

4. APPLICATION TO COSMIC STRUCTURE 

This section studies the line-correlation £(r) of simu- 
lated cosmic density fields, using both CDM and WDM, 
i.e., dark matter with a finite particle mass. 

4.1. Cosmological simulation 

We have run a suite of cosmological A^-body simula- 
tions, following the formation and evolution of LSS in a 
cubic box of comoving side- length L = 100 h~^Mpc con- 
taining 512^ ^ 1.34 • 10^ particles, fro m an initial red- 
shift o f Zinit = 199 to z = 0. Following iKomatsu et al.l 
(|2011l ) , we adopt matter and dark energy density param- 
eters of Qq = 0.273 and Qa = 0.727, a Hubble parameter 
ofh = 0.705 (defined by Hq = 100 kms-^Mpc"^), a 
primordial spectral index of ngpec = 0.95 and a normal- 
ization <J8 = 0.812. 

Initial co nditions were created using standard tech- 
niques (e.g. iPower et al.l [200 3) - a statistical realization 
of a GRF is generated in Fourier space, with variance 
given by the Hnear matter power spectrum, and the 
Zel'dovich approximation is used to compute initial par- 
ticle positions and velocities. The power spectrum for 
the CDM model is obtained by convolving the primordial 
isotropic power spectrum p{k) ex k^^p^^ with the transfer 
function appropriate for our chosen set of cosmological 
parameter s, compute d using the Boltzmann code CAMB 
fcf. iLewis e t al. 2000). Following Bode et al. (2001), the 
initial power spectra for our WDM models are obtained 
by filtering the CDM power spectrum with the transfer 
function 

TwDM(fc) = f Py^) = [1 + fc)-] , (23) 

where /i is a functi on of the WDM pa rticle mass (labeled 
'a' in eq. (A9) of iBode et all 12001 ) and u = 1.2 is a 
numerical constant. Eq. (j23j) mimics the free-streaming 
of WDM particles by preferentially suppressing high- 
frequency modes. 

All simulation s were run usin^ th e parallel TreePM 
code GADGET2 (Springel et al.""2005) with constant co- 
moving gravitational softening e = 4 kpc and individual 
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Correlation scale r/L 



Prolateness q 



Fig. 4. — (Color online) Illustration of how the function i{r) captures aspherical structure. LEFT: projection of a 3D density field 5{y) 
consisting of 50 randomly translated and rotated identical spheroids. These spheroids are characterized by the "prolateness" g, defined 
as the ratio between the pole-radius Vp and the equator-radius re, and the average radius of {rprlY^^ = 0.03L. MIDDLE: functions i{r) 
corresponding to the density fields with q=0.25 (oblate spheroids), q=l (spheres), and q=8 (prolate spheroids). In the oblate and prolate 
case, t(r) exhibits a maximum between Vp and re- RIGHT: Integral of t(r) as a function of q. 

and adaptive time steps for each particle, At = rf^feja^ 
where a is the magnitude of a particle's gravitational ac- 
celeration and T] = 0.05 determines the accuracy of the 
time integration. 

In order to study spatial correlations, the particles are 
discretized onto a regular grid of cells (here N = 400), 
as described in Appendix [Pl T he correlation functions 
are then computed via eqs. (jD4p . 

4.2. Line- correlation i{r) in a ACDM universe 

Fig. [7^ shows a projection of a simulated CDM field 
at a cosmic time of 13.7 Gyr, i.e., a redshift z = 0, while 
Fig. [Tb displays the cosmic evolution of £{r) starting at 
a cosmic time corresponding to 2:init = 199. At this early 
time, where the universe still closely resembles the ini- 
tial GRF, i{r) nearly vanishes. However, as the uni- 
verse evolves, £{r) monotonously grows for all values of r, 
hence uncovering a continuous and monotonous growth 
of phase-correlations on all scales accessible to the simu- 
lation (r < I//4 = 25/i~^Mpc). We interpret this rise of 
£{r) as a growing presence of gravity-induced aspherical 




0.001 



Filling factor / 




Fig. 6. — (Color online) Illustration of eq. (j22|) for the case of a 
3D density field composed of spheroidal functions with g = 8 (as 
in Fig. m left). The points are computational results using ro = 
0.015L (dots), ro = 0.030L (squares), and ro = 0.045L (triangles). 
Error bars represent 1-sigma shot noise uncertainties for the chosen 
numerical discretization. The solid line is the power-law of eq. (|22p 
with hi fitted to the data. 



structure, namely tidal cosmic filaments. 

We find that the late-time {z = 0) line-correlation i{r) 
is well approximated by the power-law 



^CDM(r) =0.73 



-1.15 



(24) 



0.01 

Correlation scale r/L 



Fig. 5. — (Color online) Illustration of £{r) measuring two sizes 
of filaments simultaneously. LEFT: 2D density field 6(r) com- 
posed of randomly shifted and rotated distributions e^p{—x'^/2a^ — 
y2/2cr2), where 50 distributions use (crx,cry) = (0.006L, 0.012L) 
and 10 use {ax:(Ty) = (0.02L, 0.06L). RIGHT: The corresponding 
function i(r) displays two peaks approximately centered within the 
two invervals [crx,o-y] (gray shadings). 1-sigma shot noise uncer- 
tainties of i(r) are represented by the light-red envelope. 



shown as dashed-line in Fig. [7t. A power-law behav- 
ior had to be expected for it is equivalent to i{r) be- 
ing scale- free in the sense that i{ri)/£{r2) is constant if 
ri/r2 is constant - a general feature of LSS on scales 
10 kpc < r < 10 Mpc (e.g. the weak-lensing analysis 
of the Sloan Digital Sky Survey (SDSS), see Fig. 10 of 
I Sheldon et al.l l2004f ). Interestingly, the power-law of 
eq. (|24|) is less steep than th e observed cosmic 2-PCF 
^2(^) oc r~^-^^ (jSheldon et al.l [2004). consistent with the 
interpretation that £{r) senses structural features extend- 
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Fig. 7. — (Color online) Illustration of the function £(r) in the 
case of a virtual cosmic density field, simulated as described Section 
14.11 (a) Plane-projection of the 3D field 6(r) at a simulation time of 
13.7 Gyr, i.e., at z = 0. (b) Cosmic evolution of i(r), starting from 
a Gaussian density field at z = z^nit = 199, where £{r) ~ 0. (c) 
Solid lines show the measured i(r) for four different discretizations 
of the simulation volume into A^^ cells with = 50, 100, 200, 400; 
the dashed line represents the power-law fit of eq. (j24|) . (d) i(r) for 
three different random realizations of the simulation to illustrate 
the effect of cosmic variance for a volume of (100 h~^M.pc)^ . 

ing to comparatively large scales, such as filaments. 

Fig. [7h also illustrates the numerical convergence of 
i{r) for an increasing number of grid cells (differ- 
ent solid lines). The convergence is ensured by the nor- 
malization factor in front of the integral in eq. (p!9|) as 
explained in Section [3^21 

Finally, Fig. [7]i explores the effect of cosmic variance 
on £{r). The four solid lines correspond to four different 
random realizations of the initial density field. We notice 
that cosmic variance dominantly affects the amplitude 
of the power-law, but only marginally affects its slope. 
From these examples we estimate that, for a box size L = 
100 /i~^Mpc, the normalization of 0.73 in eq. (f24|) has a 
statistical uncertainty of about 10%, while the exponent 
of —1.15 is accurate to about 3%. 

4.3. Measuring the ^temperature^ of dark matter 

What cosmological information can be extracted from 
a local measurement of i{r)7 Since i{r) only depends on 
correlations between the phases e(k), it is strictly insensi- 
tive to linear growth, defined as a uniform growth of the 
amplitudes |<5(k)|. Therefore the physics associated with 
linear growth, such as that dictating the baryon acoustic 
scale, remains invisible to £{r). However, £{r) depends 
on the physics associated with non-linear growth, namely 
on local gravitational interactions and hence on the prop- 
erties of dark matter. We therefore chose to examine 
the variations of i{r) with the dark matter particle mass 
mDM, using the WDM simulations described in Section 



Z = 100 h-^Mpc Z = 100 /z-^Mpc 

Fig. 8. — Plane-projection of two simulated 3D fields 6(r) at a 
simulation time of 13.7 Gyr, i.e., at z = 0. The two fields rely 
on identical primordial initial conditions, using CDM and WDM 
at mDM =0.1 kev, respectively. The WDM field seems smoother 
because of the suppression of short modes via eq. (|23|) . The contin- 
uous transition between the two panels reflects the periodic bound- 
ary conditions of the two boxes. 
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Fig. 9. — (Color onhne) 2-PCF ^2(r) and hne-correlation i(r) 
associated with four different dark matter particle masses. These 
examples rely on simulations with identical initial conditions. The 
data points are the galaxy-mass correlation function measured from 
weak-lensing in SDSS ( Sheldon et al...2004 ). 



14.11 An illustration of two analogous density fields with 
CDM and WDM is provided in Fig. El 

Fig.[9]shows the functions ^2(^) and £{r) for four values 
of mDM, here expressed in units of energy, e.g. 1 keV = 
1 keVc-2 ^ 1.78 • 10-^^ kg. Fig. [9] confirms that ^(r) 
only weakly depends on our array of moM, and that 
this dependance is restricted to scales r < 1 /i~^Mpc, 
in good agreement with several recent studies (e.g. 
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ISmith fc Markovid[2QTTI ). By contrast, £{r) exhibits a 
much stronger dependance on uibm , and this dependence 
extends to scales of about 5 h~^Mpc. 

The mDM-dependencies displayed in Fig. [9] can be re- 
versed to infer moM from ^2(^) and i{r). To simplify 
the notations, let /(r) be a generic placeholder for ^2(^) 
and i{r). When analyzing the variations of log^^g /(^) 
a function of moM (for r > 0.5 h~^Mpc and moM < 
1 keV), we find them to be linear in m^j^ within the 
uncertainties of cosmic variance. It follows that, at any 
r, the true particle mass can be detected against 

hyp 



a hypothetical particle mass rrijj^, e.g. against CDM 
with 

hyp\-i 



(^E^S ~ ^^^^ ^ signal-to- noise ratio proportional to 



Am 



DM 



(^DM 



«m)" 



In other words, this 



signal-to- noise ratio can be expressed as s{r)\Am-^l^\^ 
where s{r) denotes the signal-to- noise per unit of m^]^. 
We here approximate s{r) as 

/ ^ ^ I loglO /wDMQ.lkev(r) - log^p /cDM(r)| 

^^"^^ 10 keV-^ A/'(r) ' ^ ' 

where JV{r) represents the noise, defined as the standard 
deviation of log^^Q /(r) due to cosmic variance. We esti- 
mate A/'(r) as the root-mean-square of log;LO /cdm(^) over 
four random realizations of a CDM run. This numerical 
estimate is more reliable than analytical estimates based 
on the number of Fourier modes, if the latter are phase- 
correlated. The normalization factor in the denominator 
of eq. (j25j) comes from the fact that the values of ^j^J^ 
in CDM and WDM at 0.1 keV differ by 10 keV"^ 

When measuring uibm based on /(r), we consider 
Am^J^(r) to be the difference between a measurement 
of ^DM' ^ specific scale r, and its true value. 
In the Gaussian approximation, the probability dis- 
tribution (j){r) of Am^j^(r) is then proportional to 
exp [— Am^^5^(r)/2]. Note that the exponent is dimen- 
sionless, as it should be. Combining the measurements 
on the scales r G [rmin,^max] associated with indepen- 
dent Fourier modes k = 7r/r, the total probability distri- 



bution of Ar 



''DM 



becomes ( 



Hence the 



standard-deviation of a measurement of m^^l^ on scales 



r e [i 



mm? ' max. 



becomes 



cr(r, 



mm 7 ' max / 



(26) 



Numerical estimates of cr(rmin, ^max) for various inter- 
vals [rmin,^max] and for both correlation functions (^2 
and i) are given in Tab. [TJ Fig. [10] shows the probabil- 
ity distribution of Am^j^ for the case of a measurement 
based on the interval [0.5 /i~^Mpc,5 h~^Mpc]. 

The following points are worth stressing. First, the 
values in Tab. [J and probability distributions in Fig. [10] 
correspond to a perfect measurement of the density field 
(^(r), since they only account for the fundamental lim- 
itations associated with cosmic variance. They ignore 
potentially large measurement uncertainties and obser- 
vational biases. Second, the values in Tab. [1] are specific 
to our box-size L = 100 h~^Mpc. The noise JV{r) in 
eq. ([25]) scales as V~^/^, where V = is the volume of 
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Fig. 10. — Probability distribution of the difference between the 
measured inverse particle mass tu^I^ and its true value. The solid 
line represents a phase-correlation measurement via i(r), while the 
dashed line represents an amplitude-correlation measurement via 
^2(^)- Both rely on the scales 0.5 h~^Mpc < r < 5 h~^Mpc in a 
cubic volume V = (100 h-^Mpc)^ . 

the considered density field. Hence standard-deviations 
of m^l^ for any other volume V can be obtained as 
[(100 /i-^Mpc)^/V]^/V. Third, the standard-deviation 
of mDM rather than m^j^ is given by m^j^a. 

The key result from this analysis is that, in the late- 
time universe, the particle mass of dark matter moM 
is much better constrained by phase-phase correlations, 
measured via ^(r), than by amplitude-amplitude corre- 
lations, measured via ^2(^) or p{k). Furthermore, moM 
affects scales five times larger in i{r) than in ^2(^)- This 
finding is crucial, since the smaller the scale required 
to measure moM, the more the result will be entangled 
with the uncertainties of complex baryon physics, such as 
feedback fr om super novae, active black holes, and photo- 
ionization (IKim et al.ll2012[ ). 

The reason for this advantage of i{r) over ^2(^) is 
twofold. First, phase-phase correlations vanish in the 
primordial CDM/ WDM power spectrum and are strictly 
independent of linear growth. They therefore represent 
exquisite tracers of non-linear ly grown structure, which 
in turn depends substantially on moM- This supports the 
finding that £{r) is more sensitive on moM than ^2(^) for 
a given scale r. Second, £{r) is particularly sensitive to 
cosmic filaments, while ^2(^) traces spherical structures, 
such as clusters. However, cosmic filaments extend to 



[rmin,rma^]/h ^MpC 


a for ^2(^) 


a for i{r) 


[0.50,2.00] 


3.3 


0.3 


[1.00,4.00] 


> 10 


0.6 


[0.50,1.00] 


2.3 


0.3 


[1.00,2.00] 


> 10 


0.5 


[2.00,4.00] 


> 10 


1.7 



TABLE 1 

Numerical values of the standard deviations cr(rmin5 ^max) of a 
measurement of in a cosmic volume V = (100 h~^Mpc)^ . 
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larger scales than galaxy clusters and they seem to better 
preserve the primordial free-streaming scale (see Fig. [8]), 
perhaps because fil aments are secon d order effects of the 
gravitational field (jPen et al.ll2Q12f ). This might explain 
why £{r) traces moM to larger scales than ^2(^)- 

Caution is nonetheless indicated. Despite the advan- 
tage of i{r) over ^2(^)7 the case presented in Fig. [9] 
suggests that particle masses ttidm > 1 keV will re- 
main difficult to distinguish from CDM on scales larger 
than 0.5 h~^Mpc. However, particles lighter than 
1 — 2 keV seem in consistent with the Lyman-a for- 
est (|Boyarsky et al.l [2009). To probe particle masses 
^ 1 keV, £(r) will have to be me asured on smaller scales 
(eq. (3) in lSchneider et ani2012f ). not yet studied in this 
work. Those scales will still be larger than those required 
by ^2(^), but baryon physics will become important. 

5. CONCLUSION 
5.1. Three key messages 
5.1.1. i{r) - a measure of pure phase-information 

This work introduced the isotropic line-correlation 
function ^(r), defined for a density field S{r) via eq. (p!9]) . 
Unlike conventional n-PCFs, i{r) is defined exclusively 
upon the spectral phases e(k). Thus £{r) only mea- 
sures phase-phase correlations (e(k)e(q)). By contrast, 
the 2-PCF only depends on amplitude-amplitude cor- 
relations (|^(k)||J(q)|), while the n-PCFs (n > 3) de- 
pend on amplitude- amplitude correlations (|(5(k)| |(5(q)|), 

amplitude-phase correlations (|(5(k)|e(q)), and phase- 
phase correlations (e(k)e(q)). It follows that £{r) is inde- 
pendent of ^2(^) and p{k)^ up to statistical dependencies 
between the amplitudes |(5(k)| and phases e(q). 

5.1.2. i{r) - a parameter-free measure of cosmic filaments 

In Section 13.41 we have established that i{r) is, in a 
limited sense, a measure of aspherical structure, such 
as straight filaments. This measure is statistical in 
that it cannot identify the individual filaments. Hence 
typical applications of £{r) are studies of cosmic struc- 
ture rather than investigations of galaxies in particular 
environments. For the latter, sophisticated 'filament- 
finders' have been developed, i.e., algorithms able to con- 
vert 3D galaxy distrib u tions into skel etons of filaments 
(e.g. Bond et al.l 120101 : iSousbid [20TTh . On the down- 
side, the relationship between those algorithms and ro- 
bust measures of cosmic structure, suc h as co rrelation 
functions, is unknown and/or complex (i Pogosvan et al... 
I2OO9). Moreover, filament-finders always require free 
parameters, such as user-defined thresholds and scales. 
By contrast, the definition of £{r) is parameter-free and 
scale-invariant; it is even independent of the numerical 
grid, if the cell-size is much smaller than r. The rela- 
tionship between filament-finders and £{r) is therefore 
analogous to that between group-finders and ^2(^)- the 
former allow an explicit identification of individual struc- 
tural components, while the latter represent mathemat- 
ically robust measures of spatial statistics. 

5.1.3. i{r) - a thermometer for dark matter 

By calculating £{r) for simulated cosmic density fields 
(Section |4]) we demonstrated that £{r) is more sensitive 
to variations of moM than the 2-PCF ^2(^)- Therefore, 



measurements of moM are significantly more accurate 
when based on £{r) than based on ^2(^) (e.g. Fig. [TQ|) . 
Moreover, £{r) depends on moM out to scales at least 
five times larger than ^2(^)- This result is pivotal 
since smaller scales are those more affected by uncertain 
baryon physics (Kim et al. 2012) masking the footprint 
of dark matter properties. 

5.2. Prospects of using real data 

This work paves the way towards an enhanced analy- 
sis of existing an future redshift surveys, for example to 
better constrain the particle mass of dark matter moM- 

Any comparison between simulated and observed LSS 
is challenged by differences between visible matter and 
underlying dark matter. These effects are rather small 
on scales > 1 h~^Mpc. However, measurements on 
smaller scales require a precise reconstruction of the ac- 
tual dark matter density field, for example using weak- 
lensing data, and/or a modeling of visible LSS, for ex- 
ample using mock-skies based on semi-analytic modeling 
(e.g. Blaizot et al. 2005; Obreschkow et al.^2009). 

Moreover, real surveys do not come in the shape of 
a cubic box, but in a truncated survey volume with 
varying selection criteria across the volume. To deal 
with such masked data, our idealized formulation of 
the line-correlation function will need to be transcribed 
into a form applicable to a generic survey-mask. To 
do so, one might apply an approach similar to that 
of iLandv fc Szalavl () 19931 ). which essentially consists in 
comparing the correlations in the observed density field 
against those in a random field with an identical survey- 
mask. 

Finally, observational data is subject to redshift-space 
distortions (Kaiser 19871 ). leading to elongated structures 
along the line-of-sight in the reconstructed 3D density 
field ('fingers-of-God'). These spurious prolate features 
add to the line-correlation £{r). This effect will require 
additional modeling and/or an evaluation of £{r) sepa- 
rately along radial an d transverse directions, as typically 
done for the 2-PCF (jChuang Wand[20T2h . On scales 
larger than the redshift distortion scale, the radial and 
transverse parts of £{r) might also be used to constrain 
dark energy in a way analogous to the cl assical Alcock- 
Paczynski test (|Alcock fc Paczvnskilll979f ). 

5.3. Closing words 

Above all, this work demonstrates the enormous poten- 
tial of phase-information. Further investigations of this 
information may unveil a wealth of applications, extend- 
ing far beyond the case of CDM versus WDM. In light 
of future redshift surveys, the time seems ripe for phase- 
information to become a standard tool in observational 
cosmology. Ultimately, this field would tremendously 
benefit from a complete estimator of cosmic structure, 
i.e., a function J^(^(r)) that exclusively and exhaustively 
describes the information contained in a statistically ho- 
mogeneous and isotropic density field S{r). 
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APPENDIX 

A. GENERALIZED CONVOLUTION THEOREM 

We here define the Fourier transform (FT) and the corresponding inverse Fourier transform (IFT) as in iPeacockl 
(|1999 ) apart from sign of 



d{r) = lFT{d) 



V 



(27r) 



k e*''"" (5(k), 



(Ala) 
(Alb) 



where k G R^ is the wavevector. Note that the reality of 5(y) implies ^(— k) = (5*(k), where the asterisk is the complex 
conjugate. Substituting 5 for eq. ()Albp in eq. dU, we find 



Rearranging the terms, 



Sn(ri, ...,rn-i) 
Solving the integral over t gives 

Sn(ri, ...,rn-i) 



U^ 



(2^) 



3n 



n 



(A2) 



(A3) 



V 



(27r)3 



n-l 



n ^ / " \ 

i=i ^ j=i ^ 



(A4) 



where (5^ is Dirac's delta distribution in 3D. Finally, remembering that = 0, 



Sn(ri, ...,rn-i) 
which readily reduces to eq. (|8|). 



V 



n-l 



n— 1 ^ 



n-l 



(A5) 



B. EXPLICIT EXPRESSIONS FOR 6 AND ^3 

This paragraph explicits the Fourier space expressions of two isotropic correlation functions, which will be used in 
the rest of this work. First, to find the Fourier equivalent of ^2(^), we substitute S2(r) in eq. ([7]) for eq. ([12]), which 
implies 

= 7^ / cl^^ ^Dihr) P(k), (Bl) 



where w^ikr) = e^^ '^iri^r is a weighting function. A quick calculation expanded in Appendix ICl shows that 



wd{x) 



Jo(x), if I) = 2, 
sin(x)/x, if I) = 3, 



(B2) 



where Jo{x) is the 0-th order Bessel function. In eq. (|B1|) . w{kr) only depends on the integration variable k = |k|. 
The remaining D — 1 integration variables only act on P(k). Performing this integration of P(k) leads to 



6(r) 



V 



(27r) 



dk Soik) WD(kr) p(k), 

Jo 



(B3) 



where p{k) = P2{k) = |(5(k)p|jj.|^^ is the isotropic power-spectrum and Soik) denotes the surface area of the D-sphere, 
i.e., 



Soik) 



( 27rfc, if L> = 2, 
\47rfc2, ifD = 3. 



(B4) 
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Next, we consider the particular isotropic 3-PCF ^3(r) for three equidistant points on a straight line, i.e., 

6(r)^S3(r,-r)|r|^,. (B5) 
After substituting S3 for eq. (p!3|) . a derivation analogous to that in Appendix ICl then leads to 

Ur) = j^JJd''kd^qw{\k-q\r) i?(k,q). (B6) 



Note that |k — q| = y^k'^ -\- q'^ — kq cos 0^ where ^ = Z(k, q), only depends on the three coordinates (/c, 0). Thus, we 
can first integrate 5(k, q) over the remaining 2D — 3 coordinates, which leads to 



V' 



UkSnik) dqSoiq) dejoiO) 
Jo Jo 



(27r)2^yo ^^^^Vo^^^^ (B7) 

X woiV^'^ + - 2kqcos0 r) h(k,q,6), 



where h{k^q^O) = 5(k, q)n^|^^ \ci\=q z(k q)=6' isotropic bi-spectrum and joiO) is the Jacobian 



- \sin(^)/2, ifD = 3. 



C. ROTATIONAL AVERAGE OF exp(ik r) 

In two dimensions, r is expressed in polar coordinates r = |r| and ^, where 6 is the angle between k and r, such that 
k • r = /cr cos(^). Then, 

'^\r\=r = — / dOr e^^^^°^^ (CI) 
27rr Jq 

By symmetry, the real part of the integral is equal to twice the integral from to tt, and, by anti-symmetry, the 
imaginary part of the integral vanishes, 

^Iri^r = -[ dO cos{krcosO) = Jo{kr). (C2) 

Jo 

In three dimensions, r is expressed in spherical coordinates r, (/?, and ^, where is the angle between k and r, such 
that k • r = krcos{6). Then, 

'^\r\=r = I dip f dO sin l9 e^^^^°^^ = - f dO sin l9 e^^^^°^^ (C3) 

47rr Jo Jo 2 Jq 

By anti-symmetry, the imaginary part of the integral vanishes and 

'^\r\=r = - / d(9 sin i9 cos(A:rcos6>) = (C4) 
2 ./n kr 





D. NUMERICAL DISCRETIZATION 

For computational purposes, we adopt the standard numerical model: (i) the universe is described in a finite cubic 
box Q C R^ of side- length L and volume V = \Q\ = L^; (ii) this box satisfies periodic boundary conditions; (iii) 
the density perturbation field S{t) is represented on a regular Cartesian grid of cubic cells, such that the cells 
have side- lengths Ar = L/N and volumes AV = {L/N)^ . This model is valid as long as we consider correlations on 
scales larger than Ar and significantly smaller than L. The corresponding Fourier space discretization follows directly 
from the periodicity condition, which states that each mode k = (/ci, ...^ko) must satisfy kjL G 27rM Vj. Therefore, 
the Fourier cell spacing equals A/c = 2tt/L. The side- length of the Fourier box hence becomes N /S.k = 2itN/L. This 
numerical discretization is illustrated in Fig. [TT]in two dimensions {D = 2). 

We are free to choose the origins, both in direct and in Fourier space. In our convention, shown in Fig. [TTJ the sets 
of discrete vectors r and k become 

r = Ar a, where a = (ai, a^) with aj G {0, 1, {N — 1)}, (Dla) 
k = Afcb, where b = (61,..., 6^) with hj e {-fioor(N/2), -1, 0, 1, ...fioor(N/2 - 1/2)}, (Dlb) 

where fioor(x) is defined as the largest integer less than or equal to x. 
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Fig. 11. — (Color online) Discretization rules in direct (left-most panel) and Fourier (middle and right panel) space. Complex numbers 
are shown in colors as described in Fig. [2] Black dots denote our choice of the origins of the coordinate systems. 



The rules for the mapping between continuous integrals and discrete sums, both in direct and Fourier space, follow 
directly from the expressions for Ar and Ak. They read 

d^r /(r) ^ ^ Ar /(r) = E /(r)' (D2a) 

r r 

f d^fc /(k) ^ E Afc /(k) = Yl /(k). (D2b) 

k k 

wher e the value s of t he functions /(r) and /(k) on the ri ght-ha nd side are cell averages. Using the mapping rules of 
eqs. (|D2aj) and (|D2bl) . the FT and IFT of eqs. (|Haj) and (jATb]) become the discrete FT (DFT) and the inverse DFT 
(IDFT), respectively, 

r 

(5(r) = ^e^^-^ ^(k). (D3b) 

k 

By virtue of the same rules, ^2(^) in eq. (|B1[) . ^3(r) in eq. (|B6[) . and -^(r) in eq. (p^O]) become 

6(r) = ^^D(^r)P(k), (D4a) 

k 

i3{r) = J2J2'"D{\k-(i\r) B{k,ci), (D4b) 

k q 

|k|<7r/r |q|<7r/r ' V ^ 1 

These are the three functions, which we calculated in the examples, e.g. Fig. [3l 
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